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Current-Induced Pair Breaking in Magnesium Diboride 
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The transport of electrical current through a superconductor falls into three broad regimes: non- 
dissipative, dissipative but superconducting, and normal or non-superconducting. These regimes 
are demarkated by two definitions of critical current: one is the threshold current above which the 
superconductor enters a dissipative (resistive) state; the other is the thermodynamic threshold above 
which the superconductivity itself is destroyed and the superconducting order parameter A vanishes. 
The first threshold defines the conventional critical current density j c and the second defines the 
depairing (or pair-breaking) current jd- Type 11 superconductors in the mixed state have quantized 
flux vortices, which tend to move when acted upon by the Lorentz driving force of an applied 
transport current. In such a mixed state the resistance vanishes only when vortices are pinned in 
place by defects and the applied current is below the threshold j c required to overcome pinning 
and mobilize the vortices. Typically jd^>jc and a direct experimental measurement of jd over the 
entire temperature range (0 < T < T c ) is prohibited by the enormous power dissipation densities 
(p ~ I0 10 - 10 12 W/cm 3 ) needed to reach the normal state. In this work, intense pulsed signals were 
used to extend transport measurements to unprecedented power densities (p ~ 10 9 — 10 10 W/cm 3 ). 
This together with MgB2's combination of low normal-state resistivity (p„) and high transition 
temperature (T c ) have permitted a direct estimation of jd over the entire temperature range. This 
review describes our experimental investigation of current-induced depairing in MgB2 , and provides 
an introduction to the phenomenological theories of superconductivity and how the observations fit 
in their context. 

PACS numbers: 74, 74.25. Sv, 74.25. Fy, 74.25.Bt 

Keywords: Critical current, pair breaking, depairing, superconductor, superconductivity, flux, fluxon, vortex, 
magnesium diboride, mgb2 



I. INTRODUCTION 



^ ■ When a superconductor is cooled below some transition 
^temperature T c , it undergoes a phase transition leading to 
;-phe formation of a superconducting state wherein its charge 
Q^carriers correlate and condense into a coherent macroscopic 
f-*fluantum state. In the superconducting state, the system ex- 
T^j-pels magnetic flux upto magnetic field values below H c , the 
(^thermodynamic field. In type II superconductors, the flux 
**p~E>ypnkinn is partial for fields between the lower and upper 
critical fields, H c \ and H C 2 respectively. The superconduct- 
ing state is also characterized by an absence of resistivity for 
1 c urrent densities j below some critical value j c . In type II 
^superconductors, there is partial resistivity for j values be- 
tween j c and jd, where jd is the pair- breaking or depairing 
recurrent density 1 . 

| The formation of this state is governed principally by a 
• Competition between four energies: condensation, magnetic- 
^^Eield expulsion, thermal, and kinetic. The order parameter A, 
^Jhat describes the extent of condensation and the strength of 
the superconducting state, is reduced as the temperature T, 
magnetic field H, and electric current density j are increased. 
The boundary in the T-H-j phase space that separates the 
superconducting and normal states is where A vanishes, and 
the three parameters attain their critical values T C 2(H,j), 
H c %(T,j), and jd(T,H). jd sets the intrinsic upper limiting 



scale for supercurrent transport in any superconductor. 

Close to T c , jd(T) has been investigated in several 
systemsi*2i^. In a few materials^, jd(T) has been measured 
down to about 0.2T C . Here we review our investigation of 
current-induced pair- breaking in MgB2, which represents a 
complete (0 < T < T c ) investigation of jd by a direct trans- 
port method. 



1 Type I superconductors can also exhibit partial resistivity in the 
superconducting state because of fluctuations, presence of a large 
(H >H C ) current induced self field at the surface (Silsbee's rule), 
phase slippage as one approaches the phase boundary, or motion of 
magnetic domains in the intermediate state. 



II. THEORETICAL BACKGROUND 



The magnitude and temperature dependence of the pair- 
breaking current density can be adequately described by phe- 
nomenological theories of superconductivity, such as the Lon- 
don theory and the Ginzburg-Landau (GL) theory. In a mi- 
croscopic theory such as the Bardeen-Stephen-Cooper (BCS) 
theory, experimental quantities are calculated from micro- 
scopic parameters such as the strength of the effective attrac- 
tive interaction that leads to Cooper pair formation and the 
density of states at the Fermi level. Often these microscopic 
parameters are not sufficiently well known. In phenomeno- 
logical theories, connections are made between the different 
observables from contraints based on thermodynamic princi- 
ples and electrodynamical properties of the superconducting 
state. Below we give a simplified introduction to the different 
theoretical approaches for estimating the depairing current 
density and its temperature dependence. 
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A. London theory 

1. Basic framework 

The London theory^ of superconductivity provides a de- 
scription of the observed electrodynamical properties by sup- 
plementing the basic Maxwell equations by additional equa- 
tions, which contrain the possible behavior to reflect the two 
hallmarks of the superconducting state — perfect conductiv- 
ity and Meissner effect. Note that these properties hold only 
partially when mobile vortices are present. 

An ordinary metal (normal conductor) requires a driving 
electric field E to maintain a constant current against re- 
sistive losses. In the simple Drude picture, this produces 
Ohm's law behavior, j — aE, with a conductivity given by 
a = ne* 2 T/m* , where n is the number density of charge car- 
riers, and e* and m* are their effective charge and mass re- 
spectively. 

A superconductor can carry a resistanceless current and so 
an electric field is not required to carry a persistent current. 
Instead E in a perfectly conducting state causes a ballistic 
acceleration of charge so that dj/dt oc E. If the number 
density of effective charge carriers (pairs) participating in the 
supercurrent is n*, then j = n*e*v Sl where v s is the carrier 
(superfluid) velocity, and 



E = 



dj 



(1) 



This is the first London equation, which reflects the dissipa- 
tionlcss acceleration of superfluid. 

The second property that needs to be accounted for is the 
expulsion of magnetic flux by a superconductor. The mag- 
netic field is exponentially screened from the interior and 
hence follows a spatial dependence 



V 2 H = H/\ 2 , 



(2) 



where A is a screening length called the London penetration 
depth. Together with the Maxwell equation V x H = Airj /c, 
this implies the following condition between H and j: 



H 



(V x j). 



(3) 



This is the second London equation and it describes the prop- 
erty of a superconductor to exclude magnetic flux from its 
interior. The prefactors in the two London equations are re- 
lated through the Maxwell equations. Taking the curl of both 
sides of Eq.^and replacing VxEon the left side with —^^f 
and substituting for V x j on the R.H.S. with Eq. Ogives the 
expression for the penetration depth as 



A = \J to* c 2 /47m* e* 2 = \J mc 2 /4:Trn s e 2 , 



(4) 



where the pair quantities with asterisks have been replaced 
with their more common single-carrier counterparts (2n* = 
n s , e* — 2e, and m* — 2m) in the last step 2 . In the MKSA 



2 In subsequent expressions we retain the asterisk marked quantities 
in our notation to maintain generality. This allows, for example, the 
inclusion of band-structure effects on the effective mass so that a more 
exact value of m* can be used rather than simply taking the m* = 2m 
corresponding to the free-electron model. 



(SI) system, Eq.Q] becomes A = ^/10m*/47T77*e* 2 . 

A third relationship of importance in the discussion of the 
pair-breaking current concerns the thermodynamic critical 
field H c . When flux is expelled, the free energy density is 



raised by the amount H 



The critical flux expulsion 



energy (for the ideal case of a type I superconductor with a 
non-demagnetizing geometry and dimensions large compared 
to the penetration depth) corresponds to the condition 



fc = fn -fs= H 2 /8lT 



(5) 



where the L.H.S. of the equation represents the condensation 
energy — the difference in free energy densities f n — fs between 
the normal and superconducting states. 



2. Depairing current density in the London framework 

We now have all the ingredients we need to obtain the 
usual London estimate of the depairing current density. Tak- 
ing jd to represent the condition when the kinetic energy 



density equals the condensation energy, we have i 



nZm v 



*,,2 - 



\m*j 2 /(n* s e*) 2 = f = H 2 /8n. Substituting for A (Eq. HJ 
then gives the required expression for the depairing current 
density 



3d 



cHc 
4ttA' 



(6) 



In practical MKSA units, Eq. [fj] can be written as jd — 
(10 7 /47r)(i? c /A), where jd is in A/m 2 , B c is the thermody- 
namic critical flux density in Teslas and A is the penetration 
depth in meters. 

Note that this derivation assumed that n s remains un- 
changed as j approached jd- In reality n s diminishes as the 
superconductivity is destroyed and the normal phase bound- 
ary is approached. Hence Eq. [fj] will be an overestimate for 
jd- 



B. Ginzburg-Landau theory 

1. Basic framework 

There are situations where a system's quantum wavefunc- 
tion cannot be solved for by usual means because the Hamil- 
tonian is unknown or not easily approximated. The GL for- 
mulation is a clever construction that allows useful informa- 
tion and conclusions to be extracted in such a situation where 
one cannot solve the problem quantum mechanically. For de- 
scribing macroscopic properties — such as jd that we are about 
to calculate — the GL theory is in fact more amenable than 
the microscopic theory. 

The idea is to introduce a complex phenomenogical order 
parameter (pseudowavefunction) -0 = |7/>|e iV to represent the 
superconducting state. \ip{r) \ 2 is assumed to approximate the 
local density of superconducting charge carriers n s (r). The 
free energy density f s of the superconducting state is then 
expressed as a reasonable function of ip( r ) P ms other energy 
terms. A "solution" to ip(r) is now obtained by the minimiza- 
tion of free energy rather than through quantum mechan- 
ics. The unknown parameters of the theory are then solved 
in terms of measurable physical quantities thereby providing 
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contraints between the different quantities of the supercon- 
ducting state. 

Close to the phase boundary, \ij)\ 2 is small and so f s can be 
expanded keeping the lowest two orders of l^ 2 . First let us 
consider the simplest situation where there are no currents, 
gradients in \ip\, or magnetic fields present. Then we have 



2. 



fs=fn+^j\ 2 + ^\\ 



(7) 



where a and /3 are temperature dependent coefficients whose 
values are to be determined in terms of measurable parame- 
ters. The coefficients can be determined as follows. First of 
all for the solution of \ip\ 2 to be finite at the minimum free 
energy, (3 must be positive. Second, for the solution of \tp\ 2 to 
be non-zero, a must be negative. Since vanishes above 
T c , a must change its sign upon crossing T c . The minimum 
in f s occurs at 



-a/0. 



(8) 



Substituting this back in Eq.0and using the definition of H c 
(Eq.|SJ|, Eq. |S]can be recast as 



Jc 8tt 2/3 



(9) 



giving one of the connections between a and (3 and a mea- 
surable quantity (H c ). A second connection can be obtained 
by noting that n* in Eq. 0] can be replaced by \ip\ 2 , taking 
it's equilibrium value from Eq. [HI 



A = 



-0 



Aire* 



(10) 



47i#| 2 e* 2 c 

Solving Eqs. El and Ei| simultaneously gives the GL coeffi- 
cients: 

-* 2 AlTP* A 

and f3 = ff 2 A 4 (11) 



Next, the effect of fields and currents can be included in 
Eq.0by adding terms corresponding to the field energy den- 
sity and kinetic energy of the current: 



;H 2 \ 2 



~fc= <#| 







1 



a|V| 2 







w 4 



2m* 



-\ip\ 2 m*v 2 



IP 
8tt 
IP 
8tt 



(12) 



where v s is the superfluid velocity. If the amplitude is 
constant and only its phase ip varies spatially, then v s — 
/KVj£_ e_j\. \ Also in the situation of interest to us — currents 

v m cm ' 

of pair-breaking level — the gradient term is much larger than 
the A term in the operator for velocity. In this case v s w . 

It should be noted that the above derivation assumed prox- 
imity to T c only for the purpose that \ip\ be small so that f s 
could be represented as a limited power series expansion. In 
"dirty" superconductors — superconductors with a high im- 
purity scattering rate — the approximate validity of the GL 
expressions extends down to T < T c . In general, the expres- 
sions should be valid at all T as long as is small and the 
proper temperature dependent values of a and /3 are used, as 
expressed through the temperature dependencies of H c and 
A. The treatment thus far assumes that charge carriers from 
only one band contribute to the superconductivity, i.e., it is 
a SBGL (single-band Ginzburg-Landau) theory. 



current density in the SBGL framework 



We can now derive jd by finding the value of j above which 
vanishes. First we consider the case when H = 0. We will 
assume j to be uniform across the cross section. We will jus- 
tify this assumption later and look in detail at the conditions 
when the assumption is valid. Eq. 1121 then simplifies to 



2 i Ul.\A , 1 1„/.|2™*„,2 



-f c = a \^ + '-\^ + -\^ m * V 



(13) 



For zero v s , we saw earlier (Eq. [HJ that the equilibrium value 
of \ip | 2 that minimizes the free energy is |"0oo| 2 = —a/0- For 
a finite v s minimization of Eq. 1131 gives 



with the corresponding supercurrent density 



j = e*\tp\ 2 v s 







1 - 



2d 



(14) 



(15) 



The maximum possible value of this expression can now be 
identified with j d : 

(16) 



3/3 \3m* 



3V6ttA(T) 



where the GL-theory parameters were replaced by their ex- 
pressions in terms of the physical measurables H c and A 
through Eqs. ^JJ 

The approximate temperature dependence of jd can be 
obtained by inserting the generic empirical temperature 
depend encies H C (T) w ff c (Q)[l - (T/T c ) 2 } and A(T) w 
A(0)/V[1-(T/T C )4] giving 

U(T) « j d (0)[l - (T/T c ) 2 ]i[l + (T/T c ) 2 ]s (17) 

where 

j d (0) = cH c (0)/[3V6nX(0)} (18) 

is the zero-temperature depairing current density. In prac- 
tical MKSA units, Eq. [TBI can be written as j d (0) = 10 7 x 
B c (0)/[3V6ttA(G)], where j d is in A/m 2 , B c is in T and A is 
in m. 

For the dirty case, instead of Eq. El the temperature de- 
pendence of Eq. 1161 can be better approximated asi* 



3d(T)^jd(0){l-(T/T c f 



(19) 



Usually H c can't be measured reliably, but from the rela- 
tion 



<S>oH c2 
4ttA 2 



Eq. El can be recast as 



3d(0) 



c 2 $o 



216tt 3 I A 2 (0) 



(20) 



(21) 
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which in MKSA becomes jd{0) 



5.56 



10" 



\J B C 2(Q)j A 2 (0), where jd is in A/m 2 , £? C 2 is the upper critical 
flux density in Teslas and A is in meters. 

Close to T c , Eq. [T7| reduces to j d (T m T c ) m 4j d (0)[l - 
T/T c ]i. This can be inverted to give the shift in transition 
temperature T c %{j) at small currents, with the well-known 
j'2/3 proportionality: 



T c 



G) 1 


3 




Ud(P)\ 



(22) 



where T c = T C 2(0). Note that if heat removal from the sample 
is ineffective, Joule heating will give an apparent shift AT C oc 
pj 2 , which is the cube of the intrinsic ~ j 2 ^ 3 depairing shift 
near T c . Hence Joule heating is easily distinguishable from 
a pair-breaking shift. (The preceding discussion is based on 
Refs. Q and [H) 



C. GL formulation for a two-band superconductor 

We now consider the two-band Ginzburg-Landau (2BGL) 
applicable to a system such as MgB2 where two bands con- 
tribute to condensates. In this case, the condensation energy 
density can be expressed as&2iiS 



~fc = 



ailV'il 2 



"2 |^2 | 2 



02 | , ,4 



2mt 



2m* 2 





c 1 






c 1 





2 + C.C. 



"MS 



►23) 



where the first two braces correspond to the free energy con- 
tributions of the condensates of each band and the third term 
corresponds to the interband interaction energy. 

As before, for the case of no applied field and uniform cur- 
rent distribution, this simplifies to 



fc 



ail^il 2 



+ 1 - { |Vi| 2 mXt + l^lX^a} + e { ^iV>2 + c.c.} . (24) 

The phases of the two condensates are locked together since 
at equilibrium the interband free energy is minimized when 
cos(</>i — fo) = 1 or — 1 (for e < and e > respectively). 
Hence the superfluid momenta are equal, m\v\ = m^v^, and 
the superfluid velocities, and v%, will be similar to the 
extent that the effective masses are similar. 

Because of the rather large number of parameters in Ea. 1241 
it is not very meaningful to derive an expression for jd for 
direct quantitative comparison with the experimental data. 
Rather we will take the SBGL expression for j d in Eq. 
and insert the actual measured temperature dependencies of 
H c and A. We expect those empirical temperature depen- 
dencies to account for modifications due to the presence of 
two bands. As will be seen below, the experiment confirms 
this contention. Note that for current flow in the ab plane 
(which is our experimental situation) the stronger planar it 



band provides the main contribution to j s . In this case the 
behaviour should qualitatively follow SBGL with the appro- 
priately modified parameters. 



D. Depairing current from quasiparticle-energy shift 

In the previous derivations in the framework of phenomeno- 
logical theories, no account was taken of the superconducting 
gap. As a final step in obtaining theoretical estimations of j dj 
let us consider how the supercurrent is limited because of its 
effect on the superconducting gap. Unlike the GL treatments, 
this one is particularly applicable to the T -C T c regime. In 
the microscopic theory, the superfluid density does not de- 
cline continuously as j and v s are increased. Rather, n s re- 
mains roughly constant until v s reaches its depairing value 1,7 



Vd 



A 

hkp 



2h 



(25) 



which corresponds to the point when the shift, hkp.v s , in the 
quasiparticle energies (QPE) equals the gap causing the gap 
to vanish for those states. In Eq. £o is the intrinsic BCS 
coherence length. The value of £ deduced from the upper crit- 
ical field will be a lower bound on £o since scattering reduces 
the effective £ (in the dirty-limit £ w V^o^i where I is the 
mean- free pathi). j remains closely proportional to v s until 
Vd, then n s drops precipitously so that the maximum current 
is only slightly (2%) higher than the value when v s ~Vd- Thus 
to a good approximation jd — e*n*Vd- Combining this with 
Eq. |23 gives 



3d 



c 2 H 



m*c 2 A 



4ire*\ 2 hk F 27T 2 e*A 2 £ ' 



(26) 



The first R.H.S. involves the parameters such as A and Uf 
whose absolute values are not known accurately but whose 
temperature dependencies are well established, since kp is of 
course constant and the gaps have a temperature dependence 
that is well described by the standard BCS functionii*i£. 
Hence from the first R.H.S. we determine the temperature 
dependence of j d as 



3 d (T)=3 d (0) [^jvp-y 



(27) 



On the other hand, the second R.H.S. of Eq. 1261 involves 
parameters such as £o whose absolute magnitude is better 
known (from measurements of H C 2 and the relation H C 2 = 
<I>o/27r£ 2 ). Hence we use the second R.H.S. of Eq. |2S]to esti- 
mate the absolute magnitude of id(0) as 



3d(0) 



8^ 5 J A 2 (0) : 



(28) 



which in MKSA becomes jd(0) 
^AQ)/\ 2 {0). 



9.19 



io- 



E. Microscopic calculation 

Various authors have calculated jd{T) from a microscopic 
basi a 1 ! 13 ^ 4 . For arbitrary temperatures and mean free paths, 
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one must use the Gorkov equations as the starting point. 
Kupriyanov and Lukiche\ii£ have derived jd(T) from the 
Eilenberger equations, which are a simplified version of the 
Gorkov equations. This derivation is beyond the scope of the 
present review, but a nice shortened version can be found in 
reference 0. The microscopic calculation confirms the overall 
temperature dependence predicted by GL and the two nor- 
malized curves differ only slightly from each other (e.g., see 
Fig. 4 of reference 0) . 



to to be due to microscopic disorder (scattering) and not just 
some extra series resistance at grain boundaries, since H C 2 is 
also enhanced by about the same factor over its values in sin- 
gle crystals 18 ' 1 —. Also measurements of A on these films show 
no evidence of grain-boundary weak-linked behaviour, which 
would be manifested as a non-linearity in response and as an 
increase in the apparent absolute magnitude of A, which were 
not observed^. 



F. Comparison between different theoretical estimates 

We will now compare the different theoretical estimates of 
jd(0) obtained by the different approaches and calculate their 
relative ratios from Eqs. H3 IS] an d ESI The London es- 
timate, as discussed earlier, should be an overestimate. Its 
value is 1.84 times higher than the GL estimate. Similarly, 
the estimate from QPE shift turns out to be 1.65 times higher 
than GL. Since the QPE shift calculation is based on a sim- 
ple single isotropic gap, the actual jd(0) will be different and 
somewhat lower. For the very dirty (/ <C £o) an d very clean 
{I ^ £o) limits, the actual jd(0) is expected to be 0.67 and 
0.82 times the London estimate respectively 1 . Hence these 
estimates should be correct at most to a factor-of-two accu- 
racy. Errors in the values of the parameters will further add 
to the inaccuracy of the calculated value. So when comparing 
with the experiment, an agreement to within half an order of 
magnitude is about the best that can be expected. 

III. EXPERIMENTAL METHODS 
A. Samples 

The samples are 400nm-thick c-axis oriented films of MgB2 
on sapphire. In this paper we show data on four bridges, la- 
belled S, M, N, and L with lateral dimensions 2.8 x 33, 3.0 
x 61, 3.0 x 60, and 9.7 x 172 /jm 2 respectively. The lateral 
dimensions are uncertain by ±0.7/im and the mean thick- 
ness by ±50 run. TEM (transmission electron microscopy) 
shows variations in thickness of about 8t/t ~ 10%. The films 
were fabricated using a two-step method whose details are 
described elsewher o 16 ! 17 . An amorphous boron film was de- 
posited on a (1102) AI2O3 substrate at room temperature 
by pulsed-laser ablation. The boron film was then put into 
a Nb tube with high-purity Mg metal (99.9%) and the Nb 
tube was then sealed using an arc furnace in an argon atmo- 
sphere. Finally, the heat treatment was carried out at 900° 
C for 30 min. in an evacuated quartz ampoule sealed un- 
der high vacuum. X-ray diffraction indicates a highly c-axis- 
oriented crystal structure normal to the substrate with un- 
detectable (< 0.1%) impurity phases. Magnetization M(T) 
curves have a A-limited transition width of 1.5K (T c -spread 
< 0.2K). However, the aforementioned variations in thickness 
produce a broadening of R(T) with increasing j. The films 
were photolithographically patterned down to narrow bridges 
and in the case of one sample (N) the lead areas were further 
delineated by mechanical scribing. 

The samples have a normal-state resistivity at room tem- 
perature p„(300K) sa 14^f2-cm, which is about 7 times that 
found in clean single crystals 18 . The enhancement appears 



B. Cryogenics 

Two sets of apparatus were used for the measure- 
ments. Samples S, M, and L were measured in an Oxford 
Instruments^- 1 liquid-helium based vapour-flow cryostat with 
a 16 Tesla superconducting magnet. Sample N was measured 
in a Cryomech^ pulsed-tube closed-cycle refrigerator and 
a water-cooled copper magnet. In the latter case, the "zero 
field" could be made very small (< 0.7 G) by shielding the 
sample region using mu-metal. Lakeshore^ diode and cer- 
nox sensors were employed for the thermometry. 

C. Electrical measurements 

Low-current resistance measurements (/ < 50/xA) were 
made by the standard four-probe DC technique. Higher cur- 
rent measurements were carried out using a pulsed technique 
with a four-probe configuration. A pulsed signal source (ca- 
pable of both constant-current and constant-voltage modes) 
drives the signal through the sample and a series standard re- 
sistor R-gtd- The signals across the sample and Rg^d ar e then 
detected with a digital-storage oscilloscope. The temporal 
reproducibility of the system is ~1 ns as can be evidenced 
from Fig. 2] which shows two sets of measurements of I(t) 
and V(t). The relative delay between rising edges of I it) and 
V(t) corresponds approximately to the extra length of wires 
divided by the speed of light. 



0.08- 
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FIG. 1: Rising edges of voltage and current (voltage across R s td) 
across an MgB2 bridge (sample N). The sample is in the nor- 
mal state at room temperature. The graph contains two sets of 
measurements of I(t) and V(t) showing overlap and temporal re- 
producibility to within about 1 ns. 

Pulse durations range 0.1-4 /is with a duty cycle of about 
1 ppm. About 100 pulses are averaged to reduce the noise. A 
typical pair of current and voltage pulses is shown in Fig. [21 
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FIG. 2: Pulse waveforms at high dissipation levels (j = 9.7 
MA / cm 2 , E = 83 V/cm, and p = jE = 803 MW/cm 3 on the 
plateaus). The measurement was done on an MgP)2 bridge (sam- 
ple S) just above T c (at 42K). The resistance rises to (90% of) its 
final value in about 50 ns from the (10%) onset of I. 



The values of I and V used in subsequent analyses correspond 
to the flat plateaus of the pulse waveforms. As these are time 
independent, the measurement corresponds effectively to a 
DC measurement albeit over a shorter-than-usual duration. 
The computed R(t) = V(t)/I(t), shown in Fig. is seen to 
have a 50 ns rise time. The total duration of the pulse is 
not relevant but the time at which the voltage or current is 
measured since it turns on. Since the signals saturate to well 
defined values by — 100 ns, this is the effective dissipation time 
r that is relevant to the thermal calculations below. Some 
additional information on the electrical-measurement setup 
can be found in a previous review article 21 and other recent 



D. Heat conduction during pulsed measurements 

We now look in more detail at the thermal processes in- 
volved in the removal of heat from the sampl o 24 ' 2 ^ . Heat 
generated in the film diffuses toward the interface with the 
substrate essentially instantaneously. On the time scale of 
nanoseconds, phonons transfer heat across the interface be- 
tween the film and substrate. Heat then diffuses within the 
substrate in a matter of microseconds and finally into the 
heat sink in milliseconds. Thus for pulses of microsecond or 
shorter duration, the heat doesn't leave the substrate and fol- 
low duty cycles the thermometer will not register a temper- 
ature rise. The sample temperature rise is then composed of 
a temperature gradient within the film, an abrupt temper- 
ature difference across the boundary, and a gradient within 
the substrate. Each can be expressed as an additive compo- 
nent of the thermal resistance Rth, which we will define as 
the temperature rise AT per unit power density p. 

The temperature variation within the film given by 
ATi(z) ~ (p/nf)(tz — z 2 /2) where (t— z) is the distance of a 
point within the film from the interface, t is the film thickness 
and Kf is the thermal conductivity of the film. Expressed as 
a thermal resistance this becomes 



Rthi(z) 



1 



tz 



(29) 



For MgB 2 , K f = 0.1 and 0.4 W/K.cm at 10 K and 40 K 
respectively 2 ^**^. This yields an average R t h in the midplane 
of the film (i.e., at z = t/2) of 6 nK.cm 3 /W at 10K and 1.5 
nK.cm 3 /W near T c . 

Phonon mismatch at the interface between film and sub- 
strate produces an abrupt temperature drop AT2 = ptRbd, 
where Rm is the thermal boundary resistance at the film- 
substrate interface. This can be expressed as a second com- 
ponent to Rth- 



Rth2 = tRbd- 



(30) 



Most metallic films on sapphire exhibit the approximate em- 
pirical rule R bd T 3 a 20 cm 2 .K 4 /W 28 . Thus R th2 - 12 
nK.cm 3 /W at 40 K (near T c ) and -800 nK.cm 3 /W at 10 
K. (Ceramic on ceramic boundaries, such a YiBagCugOy on 
LaA103 seem to fair better in this respect.) 

Finally as the heat pulse propagates through the substrate, 
there is a third component to Rth , which for the long rectan- 
gular strip geometry is given by 2 ^: 



Rth3 = 



2.26trw 



2(D s t) 1 / 2 [4(D s t) 1 ^ + w]c s 



(31) 



and c p are the 



where D s — n s /c s is the diffusion constant, k s 
thermal conductivity and specific heat of the substrate ma- 
terial, w is the width of the bridge, and r is the duration for 
which the power is applied. For sapphire, c s — 3.2 J/K.cm 329 
and K s» 10 and 100 W/cm.K at 10 K and 40 K respectively^. 
Thus D s — k s /c s — 3 and 30 cm 2 /s at 10 and 40 K respec- 
tively. This gives a thermal diffusion lengths \J Dt within 
the substrate of 6-17 (tzm for T=10-40 K, for r = 100 ns. 
Since our bridge widths (—2-10 /an) are small or comparable 
to these diffusion lengths, the denominator — 8Dtc s cx k s is 
mainly dependent on the conductivity of the substrate and 
is relatively independent of its specific heat. From Eq. |3"T1 
the last component of the thermal resistance then has an es- 
timated magnitude of R t h3 ~ 0.8-0.1 nK.cm 3 /W in the T 
range 10-40 K respectively for a 3 /zm wide bridge. 

From the above estimated components, we can obtain the 
total thermal resistance 



R 



th 



R 



thl 



R, 



th2 



Rt 



//.'! 



(32) 



as R th = 1.5 + 10 + 0.1 « 12 nK.cm 3 /W near T c and 0.8 
/iK.cm 3 /W at 10 K, where the dominating term comes from 
the boundary resistance. 



E. Adiabatic heating 

The worst-case scenario for sample heating is when the 
timescales and conditions are such that none of the generated 
heat escapes. In this case the energy density pr dumped by 
the pulse will go entirely into raising the internal energy U. 
The temperature rise can then be expressed in terms of the 
film's specific heat as 



pr = SU 



( c(T)dT « cAT, 

JTo 



(33) 



where the last approximation applies for small temperature 
shifts, c = c(Tq) is the specific heat at the nominal bath 
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FIG. 3: Pulse waveforms showing a time-dependent temperature rise at extreme dissipation levels (p > 10 10 W/cm 3 ). The measurement 
was made on a 100 nm thick YiBa 2 Cu 3 07 film on LaAlOj^ 3 . a t T=50 K and B=10 T. The left, middle, and right panels show the 
current density, electric field, and power-dissipation density respectively. 



temperature To, and T' is the final raised temperature, i.e., 
AT = T' — T). The effective thermal resistance in this case 
is 

R thA = r/c(T ). (34) 

Close to T c , c(40K) w 0.044 J/K.cm 3 in MgB^. For times 
between 100ns and 1/Lts after current turn on, the corre- 
sponding adiabatic Rth then has values in the range ~2- 
20 /iK.cm 3 /W. Note that the RthA values are about three 
orders of magnitude higher than when there is heat conduc- 
tion (Eq. I32[l but place an absolute upper limit on the actual 
Rth- Because RthA is directly proportional to r, the amount 
of heating can be made almost arbitrarily small by reduc- 
ing the duration of the signal (for very short durations there 
will be addional corrections when the electrons can no longer 
equilibrate with the phonons). The above estimates were for 
T~T C , at 5K the corresponding RthA values will be one or 
two orders of magnitude higher. 



F. Heat generation at contacts 



G. Pulse waveform distortion due to temperature rise 

The third component of the thermal resistance Rth?, arising 
from the flow of heat into the substrate is time dependent 
(Eq. 1311) . As a result intense heating can cause the resis- 
tance to rise with time and hence the voltage pulse becomes 
slopedSi. Fig. [3| shows an example of this situation, where 
the measurements were done on YiBa2Cu307 films. Even in 
the case of a system with a relatively T independent p n such 
as the MgB2 films, in the superconducting state the flux-flow 
resistivity is T dependent, especially in the regime near T c . 

For very short pulses and/or low enough power densities, 
the temperature rise in the substrate AT = pRth3 becomes 
negligible and the main bottle neck for heat conduction be- 
comes the thermal boundary resistance. If these conditions 
are satisfied the voltage pulse should appear flat. This was 
verified in all of our measurements and data was only taken 
when the pulses were undistorted as in Fig. |2 Note that the 
power-dissipation densities for the pulses in Fig. |3| are much 
higher than the levels used in the present work on MgB2. 



As seen above, pulsed measurements greatly reduce the 
effective thermal resistance and accompanying temperature 
rise in the sample. In addition they can essentially elimi- 
nate problems associated with heat produced at lead con- 
tacts. This is because the contacts are typically located sev- 
eral mm away from the active part of the sample (the bridge) 
involved in the four-probe measurement. This distance is 
typically large compared to the thermal diffusion distance 
y/Dr. The specific heat at 10 K and 40 K is 1,400 and 50,000 
J/K.m 3 respectively^ 3 ^. Together with the values of 10 and 
40 W/m.K for thermal conductivity^*^ gives diffusion con- 
stants of 0.007 and 0.0008 m 2 /s at the two temperatures. 
Thus for dissipation durations 0.1-1 [is and temperatures 10- 
40 K, we get diffusion distances in the 9-84 /im range, which 
are about a hundred times shorter than the distance to the 
bridge. So heat generated at the contacts does not interfere 
with the measurement. 



H. Influence of thermal environment 

In measurements involving a continuous DC signal or long 
time scales, substantial overall heating can be alleviated by 
immersing the sample directly in a liquid cryogen. This is es- 
pecially noticeable when the cryogen is in a superfluid stated. 
In our measurements the timescales are such that even the 
substrate does not experience a sigificant rise in T. In this 
case we would expect that additional heat removal from the 
exposed surface of the film to a cryogen to not make a big dif- 
ference. This is indeed the case as can be seen in Fig.0]where 
IV curves measured in different thermal environments show 
a jump to the normal state at a value of Id that is roughly 
independent of the thermal environment. 
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FIG. 4: IV curves for sample L at the lowest temperatures with 
the exposed surface of the film in contact with different thermal 
environments as indicated. The influence of thermal environment 
on Joule heating is seen to be minimal as there is no major system- 
atic change in the threshold current at which the jump in voltage 
occurs. 



occurs mainly because of a drop in H C 2 with T). This will 
lead to an insignificant rise in pf even for a constant current. 
If the measurement is made with a constant-voltage source, 
Pf will actually drop with R (since then pf oc V 2 /R instead 
of pf oc I 2 R). Our apparatus, when set to constant-voltage 
mode, has a source impedence as low as 0.3 Ohms which is 
about an order of magnitude below R n . Thus the power dis- 
sipated doesn't rise rapidly with T. In the meanwhile the 
thermal properties that aid heat removal (conductivities and 
specific heat) become rapidly more effective as T increases, 
working strongly to prevent a runaway. Thus it is possible 
to make a quantitative evaluation of sample heating and the 
risk of thermal runaway. 

A more detailed discussion on the subject of hot spots and 
thermal propagation in superconducting microbridges can be 
found in the review article by Gurevich and Mints^ 3 - and ref- 
erences therein. 



K. Depinning current versus depairing current 



I. Electron-phonon disequilibrium 

Beside macroscopic heating of the sample, other exotic sit- 
uations can arise where the electronic temperature is signif- 
icantly raised with respect to the phonons, which remain 
in equilibrium with the bath 23 i 34 i 35 i 36 i 37 i 38 i 39 i 40 i 41 i 42 i 43 i 44 i 45 . 
The characteristic time for the relaxation of energy between 
the electrons and phonons is typically under 1CP 8 s. This 
effect happens mainly when there is vigorous vortex motion 
and is not a main concern in the present work. 



As long as gross sample heating or thermal runaway can be 
eliminated as discussed above, there is no confusion between 
the conventional (depinning) j c and jd- When jd is exceeded, 
the resistance reaches the normal-state value R n . When j c is 
exceeded, the resistance approaches (if pinning is well over- 
come) the free- flux- flow value Rf ~ R n B/H C 2. The latter is 
lower than R n by three orders of magnitude. Hence the crit- 
ical current measured is not just the conventional depinning 
threshold but a good estimate of jd- 

L. Uniformity of current flow 



J. Temperature rise and thermal runaway 

Except in the case when the sample has perfect electri- 
cal conductivity, in the dissipative state the measurement is 
necessarily accompanied by some rise in sample temperature. 
This leads to the risk of thermal runway, where a small tem- 
perature rise leads to increased dissipation, leading to a fur- 
ther temperature rise and so on. The sample can end up at 
a temperature much higher than the bath temperature and 
possibly in the normal state. If this scenario is fulfilled, the 
quantity actually measured will be the conventional j c , de- 
mar kating the onset of dissipation, rather than the true jd- 
The match needed to light a thermal runaway is flux-flow dis- 
sipation whose power density is given by pf < j 2 p„B / H C 2- A 
small localized hot spot will not sustain a runaway since the 
current can simply avoid the hot spot and flow through more 
conductive regions. Also heat produced at one point does 
not remain localized on the timescale of the measurement 
but spreads uniformly across the width of the film because 
the thermal diffusion length, calculated earlier, is larger than 
the bridge width. Thus the the relevant parameter is the 
average dissipation across the sample volume. 

At 10 K, this pf is about 1.2 W/cm 3 taking an average 
self field of 100 G (all measurements here were made in zero 
applied magnetic field with the sample chamber shielded by 
mu- metal in some cases). From the Rth calculated earlier, 
this results in AT sa 1.1 K, which increases the sample re- 
sistance only marginally (p n of MgB2 is fairly flat; the rise 



A superconducting wire with dimensions small compared 
to £ and A, automatically has a uniform distribution of j and 
|A| across the cross section. Close to T c , £ and A diverge, so 
this condition applies to the T C 2 shifts measured from resistive 
transitions at fixed currents (Eg. 1221) . 

At lower temperatures and for sample cross sections that 
are not small, the question of current-flow uniformity needs 
to be further examined. In a type-I superconductor in zero 
applied field, the self flux is completely expelled and the cur- 
rent flows without resistance along the periphery of the cross 
section up to some threshold value. Resistance first appears 
when the self field at the surface exceeds H c (Silsbee's rule). 
Beyond this the sample enters an intermediate state where 
normal regions coexist with superconducting ones and cur- 
rent flows through both so that the macroscopic resistance 
is a fraction of the normal-state value. In this intermediate 
state, the current is distributed more homogeneously over the 
sample volume. 

In thin strips or bridges of type-II superconductors this 
process goes a step further. First, instead of H c you have 
a much lower H c \ that defines the threshold of flux entry 
(possibly modified by surface barriers except for thin film 
bridges). Likhare^^ has shown that the minimum sample 
width for the nucleation of vortices is 4.4£(T). Second, be- 
cause of the highly aspected geometry and consequent large 
demagnetization, the effective threshold is practically negli- 
gible. Thus there are always flux filled regions within bridges 
wider than the Likharev threshold. Also for the type-II case, 



9 



^10 

E 
o 

G 

a 1 



□ MY (edge) 
• MY (uniform) ch~ 

x ly(-) D a D „y&r 



□ 
□ 



x i 

X 



p~1%p n 



1 

j (MA/cm 2 ) 

FIG. 5: Non-linear mixed-state response of two 90-nm thick 
YiBa2Cu307 microbridges. Samples MY and LY are 8 /im and 
16 pm wide respectively. The lower two sets of symbols, showing 
p and j calculated assuming uniform current flow, are in excel- 
lent agreement. The open squares at the top, showing p and j 
(for MY relative to LY) based on an assumption of non-uniform 
current flow along edges, are in conspicuous disagreement. The 
dashed line corresponds to 1% of the normal-state resistivity. The 
measurements were made at T = 50 K and B — 13. 8T (~ 20% of 
H c2 ); A ~ 200 nm. 



IV. RESULTS AND ANALYSIS 



A. Resistive transitions at fixed currents 
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vortices are present that necessarily move since we are con- 
cerned with current densities well beyond the depinning j c . 
Hence the current flow will become homogeneous under these 
conditions. In general, as the current grows beyond j c and 
the system becomes highly resistive, the current flow becomes 
macroscopically uniform, as in a normal conductor, due to 
the principle of minimum entropy production. This has been 
investigated by us in our previous work on YiBa2Cu30y 
bridges of different widths 22 . Fig. [5] shows p{j) curves on 
two YiBa2Cu307 bridges measured in the highly non-linear 
mixed-state regime. Note that the resistivity changes by al- 
most two order of magnitude and extends below 1% of p n . 
The values of p and j calculated based on uniform current 
flow (i.e., dividing / by the nominal cross section to obtain 
j, etc.) are in excellent agreement. On the other hand if 
one assumes that the current flows non-uniformly along the 
edges, the nominal cross section will overestimate the actual 
area that contains the current. In this case, the relative error 
between samples MY and LY will be a factor of two in p and 
a factor of half in j (the sample widths differ by a factor of 
half) . The edge flow assumption is checked by comparing the 
squares to the crosses. Clearly the assumption fails. We con- 
sider this to be definitive proof that the current does not just 
flow along the edges but occupies the whole sample cross sec- 
tion uniformly once the transport becomes dissipative. Other 
authors^i have also reached the same conclusion that vortex 
motion homogenizes current flow. 

It will be seen in the results below that the values of jd(Q) in 
MgB2 obtained from resistive-transition shifts near T c agree 
with the values obtained at low temperatures from IV jumps. 
Also the values obtained for samples with different cross sec- 
tions are all in mutual agreement. If the current flow did not 
fill the cross section, the macroscopically calculated value of 
jd would be higher for a narrower sample. All of these ob- 
servations confirm that current flow becomes homogeneous 
under dissipative conditions. 



FIG. 6: Resistive transitions of MgE>2 bridges at different currents 
(values correspond to curves from left to right.). Panels (a) and 
(b) show two windows of the same data. The inset in (a) shows 
the sample geometry and configuration of leads. Panels (b), (c), 
and (d) show the central main portions of the transitions for three 
different sized samples. The rightmost curves at 1=1.4 pA were 
measured with a continuous DC current; the rest used pulsed sig- 
nals. The dashed lines represents R = R n /2 for each sample. 
Reprinted with permission from M. N. Kunchur, S.-I. Lee, and 
W. N. Kang Phys Rev. B 68, 064516 (2003). (c) (2003) The 
American Physical Society. 

Fig. EJa) shows the resistive transitions at different elec- 
tric currents I for sample M. The inset shows the sample 
geometry. The horizontal sections of the current leads add 
a small (~ 15 %) series resistance to the actual resistance 
of the bridge. Because j in these wide regions is negligi- 
ble, this resistance freezes out at the nominal unshifted T c , 
making the onset seem to not shift. Similarly, the lower 
foot of the transition will have a flux-motion contribution 
Rf ~ R n B/H c2 < 5% R n from the self field. The central 
two-thirds portion of the transitions (magnified in panel (b)) 
circumvents these errors, displaying relatively parallel shifts 
due to pair breaking. 

Variations in film thickness cause the transitions to 
broaden slightly with increasing j despite phase purity For 
a simple model with series thickness variations St, the func- 
tional shape of the current-dependent broadened transition 
is given by R(j,T,5t) = R n {log(j /Aj d (0)) - 1.51og(l - 
T/T c )}/ log(l + St/t), where <ft=thickness variation. At the 
Rn/2 criterion (shown by the dashed line) the actual T c shifts 
correspond exactly to shifts for a sample with the same mean 
thickness t but with St = 0. 

Panels (c) and (d) show similar sets of curves for two more 
samples. It should be noted that in addition to transition 
broadening due to thickness variations, the R(T) transition 
may have some intrinsic width as a function of j. An analo- 
gous situation arises for R(T) in a magnetic field, where there 
is an intrinsic broadening as B is increased^. Unfortunately 
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there is no theoretical work on the R(T) transition shape at 
high j. Nevertheless, we expect the midpoint criterion for the 
shifted T c to provide a factor-of-two estimate of jd- Fig- da) 
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FIG. 7: Shifted transition temperatures at different currents. The 
two panels compare the same T C 2(j) data versus 7 2/3 and I 2 , show- 
ing adherance to the I 2 ^ 3 law for pair breaking rather than the I 2 
law for Joule heating. The linear fit (solid line) to the I 2 ^ 3 plot 
gives I d (0) = 257 mA (see Eq-HH- 

shows the midpoint T^s and their corresponding currents 
(ranging from 10~ 6 to 1(T 2 A) plotted as I 2/3 (expected for 
pair- breaking) and as I 2 (expected for Joule heating). The 
shifts are closely proportional to I 2 l z rather than to I 2 , show- 
ing that heating is not appreciable (the plots for the other 
samples look similar). In fact from our earlier calculation of 
the thermal resistance, we can estimate the expected temper- 
ature rise. For the average j m 0.4 MA/cm 2 , and transition- 
midpoint p w 3/ifl-cm, we get p w 0.5 MW/cm 3 . This heats 
up the sample by AT = pRth ~ 6 mK. This is much smaller 
than the size of the symbols on the plot and two orders of 
magnitude smaller than the observed shifts. Therefore heat- 
ing is definitely negligible for the data shown in Fig. 

The slope dl 2 / 3 /dT c (j) together with Eg . |2*2*I gives a zero- 
temperature depairing current value of 257 mA (if the T c 
criterion is taken at 30% and 70% of R n , the corresponding 
Id values are 196 mA and 299 mA respectively). Dividing this 
by the cross-sectional area gives jd(0). For samples S, M, N, 
and L the respective values of jd(0) are 2.2 x 10 7 , 2.1 x 10 7 , 
2.0 x 10 7 , and 1.8 x 10 7 A/cm 2 . The four values are consistent 
within the uncertainties in the sample dimensions, implying 
a cross-sectionally uniform current density. As discussed ear- 
lier, this is especially expected close to the T C 2(j) boundary 
where A and £ diverge and the superconductor becomes highly 
dissipative due to flux motion and fluctuations. 
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FIG. 8: IV curves for sample N at fixed temperatures (listed 
for curves going from right to left). Lines are drawn to guide the 
eye. Curves at intermediate temperatures were omitted for clarity. 
Beyond some threshold, the voltage jumps abruptly to the normal 
state. The slanted portions have the slope V/I = R„ and their 
intercepts are zero. 



So the measured threshold can be associated with Id(T) at 
the nominal bath temperature. 

For T > 35 K the transition is gradual whereas at the lower 
temperatures it is rather abrupt 49 . This is partly because of 
film thickness variation as discussed earlier and partly be- 
cause a type II superconducting phase transition changes 
from second order to first order at lower temperatures in the 
presence of a current. 

The measurements on sample N were done with a signal 
source with a high source impedance thereby eliminating the 
"s" shape seen earlier in reference^, where the external cir- 
cuit had a source impedance comparable to that of the sam- 
ple. 



C. Zero-temperature values of jd 

From such IV characteristics measured at the lowest tem- 
perature the current required to drive the sample normal pro- 
vides a direct estimate of jd{T w 0). These values of jd(0) for 
the four samples are shown in Table 1 (top row of numbers) . 
These are seen to be consistent with the values (bottom row 



Method: 


Sample 


S 

1.9 x 10 7 
2.2 x 10 7 


M 

2.0 x 10 7 

2.1 x 10 7 


N 

2.3 x 10 7 
2.0 x 10 7 


L 

1.7 x 10 7 

1.8 x 10 7 


IV at T < T c 

R(T) shift at T ~ T c 



B. Current-voltage curves at fixed temperatures 

Fig.[S]shows current-voltage (IV) characteristics at various 
fixed temperatures for sample N (results for samples S, M, 
and L are similar). As / is increased, V remains close to zero 
until some threshold value. Above this it switches abruptly 
to the Ohmic behaviour V — IR n . This threshold is a lower 
bound for Id', however, as per our earlier discussion regarding 
thermal runaway, the temperature error is of the order of 1 K 
at low temperatures and negligible at higher temperatures. 



Theory: 




GL 


4 x 10 7 


QPE shift 


7 x 10 7 



TABLE I: Zero-temperature values of jd in A/cm 2 

of numbers in Table 1) obtained from shifts in the resistive 
transitions near T c (Fig.|7|and Eg. I22fl . It may be reiterated 
that the measurement does not reflect a depinning j c ; with- 
out significant pairbreaking or runaway heating, the motion 
of the minuscule self flux (B <C-ff C 2) will produce Rf <C R n - 
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The observed R ss R n is reached only when the supercon- 
ductivity is completely destroyed and the system has become 
normal. 

The average value for all samples by both methods is 
j d (0) « 2± 0.7 x 10 7 A/cm 2 . 

Our experimental estimate of jd(0) can be compared with 
theoretically calculated values from Eq. |^ and Eq. |2S] The 
H C 2 that enters these equations is the one that reflects the 
clean-limit BCS coherence length in the ab crystalline plane 
(H || c axis). An actual measured H c2 reflects the reduced 
coherence length due to scattering (£ Vfo^)- Hence jd cal- 
culated from an empirically measured H C 2 is going to be an 
overestimate. We take H C 2(0) ~ 3 T obtained from crystals 
by Sologubenko et as an upper bound on the clean-limit 
H C 2- Because of sample to sample variation and different 
amounts of impurity scattering, the uncertainty and range in 
measured values of H C 2, span a factor of five^S. For A we take 
the value of 150 nm found using an AC induction technique 
on the same type of film samples by Kim et al^S. The exper- 
imental uncertainty and range in A is about ±40%£i, which 
is a tighter error bar than that of H C 2- With these values 
of H C 2 and A we obtain jd(0) based on the GL calculation 
(Eq. I21|) and the quasiparticle-energy shift (Eq. I28fl . These 
calculated values are also shown in Table 1 and are indeed 
higher than the measured jd (by about a factor of 2-3), but 
are of the same order of magnitude, which is about the best 
agreement that can be expected given the uncertainties in 
the parameters. 

Incidentally, from the measured valued of 150 nm and 
Eq.0J the superfluid density turns out to be n s (0) « 1.3xl0 21 
cm~ 3 and the measured^ H C 2(0) ~ 3 T and Eq. 1201 give 
H c « 0.18 T. 
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FIG. 9: Pair-breaking currents from IV characteristics at fixed 
temperatures (jumps in Fig.|HJl. The dashed line represents the GL 
Id(T) function (Fq. ll7H and the solid line represents the theoretical 
Id{T) resulting from quasiparticle-energy shift (Eq. I27H . The end 
points ld (0) and T C 2(0) were fixed by the data; no other adjustable 
parameters are involved in the interpolation. 



D. Temperature dependence of Id 

We now take the Id values determined by the IV jumps 
of Fig. [5] and plot them versus temperature. Fig. [5] shows 
such plots for two samples. The lines show the theoreti- 
cal temperature dependencies expected from the GL treat- 
ment (Eq.Uni and from QPE shifts fEa.l27|). For H C (T) we 
take the generic empirical temperature dependence H C (T) sa 
i? c (0)[l - (T/T c ) 2 }, for A(T) we take the BCS function, and 
for A(T) we take the empirical temperature dependence mea- 
sured for MgB2 by Kim et al. 3 . The ends of the curves 
[0,7(2(0)] and [T C 2(0),0] are fixed. Other than that there are 
no adjustable parameters. As can be seen the data follow the 
general trend of the theoretical curves. 

V. CONCLUSIONS 

In conclusion, we have studied current induced pair- 
breaking in magnesium diboride over the entire temperature 



3 For our purpose and temperature range of interest (0 < T < 0.95T C ), 
Kim et al.'s empirical temperature dependencaSfi of A can be suf- 
ficiently well approximated by the simple function A(T) m A(0)[1 — 
(T/Tc) 2 ' 2 ]" 1 / 2 forT < 0.8T C and A(T) ss 2.74A(0)[1 - (T/T c ) 1 ' 2 ]" 1 /2 
for T > 0.8T C . Similarly the BCS temperature dependence of the gap 
can be approximated by A(T) « A(0)[1 - tan(0.67396 x (T/T c ) 3 7 )]. 



range for in-plane current transport. The measured Id{T) 
function is consistent with the expected theoretical temper- 
ature dependencies and conforms exactly to the AT C cx j 2 / 3, 
behavior predicted near T c . jd{0) obtained from the value 
of current required to drive the sample normal at T — > 0, 
agrees with the ^(0) deduced from the AT C cx j 2 / 3 behav- 
ior close to T c . The average value from our measurement 
is jd(Q) ~ 2 ± 0.7 x 10 7 A/cm 2 , which is about a factor of 
three lower than the value ~ 6 x 10 7 A/cm 2 calculated from 
published parameters, but the two are consistent within the 
uncertainities of the various parameters. 

From a technological standpoint, the depairing current 
density of MgB2 is about an order of magnitude lower than 
the high-T c cuprates^. The good news is that flux motion in 
films is so quenched^ that the depinning j c at modest fields 
appears to be over 25% the magnitude of j^^^, whereas 
for the cuprates, j c and j d can be separated by two or three 
orders of magnitude^. 

The tremendous experimental difficulties against measur- 
ing j d (0) until now, can be appreciated when one sees that for 
YiBa2Cu307 (where j d (T w T c ) was measured^) the power 
density would be 2 ^ pj 2 - 10~ 4 (10 8 ) 2 - 10 12 W/cm 3 !— 
hopelessly beyond our pulsed technique's limit of ~10 10 
W/cm 3 . Low-T c materials like Nb and Pb also have pro- 
hibitive pj 2 values. MgB2's parameters (pj 2 ~ 10~ 5 (10 7 ) 2 ~ 
10 9 W/cm 3 ) brought jd(0) within experimental reach. 
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